Decoherence and the rate of entropy production in chaotic quantum systems 
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We show that for an open quantum system which is classically chaotic (a quartic double well 
with harmonic driving coupled to a sea of harmonic oscillators) the rate of entropy production 
has, as a function of time, two relevant regimes: For short times it is proportional to the 
diffusion coefficient (hxed by the system-environment coupling strength). For longer times 
(but before equilibration) there is a regime where the entropy production rate is fixed by the 
Lyapunov exponent. The nature of the transition time between both regimes is investigated. 
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Environment induced decoherence has been identified 
in recent years as one of the main ingredients in the tran- 
sition from quantum to classical behavior. Classicality 
emerges as a consequence of the coupling of quantum 
systems to an environment which, in effect, dynamically 
enforces super-selection rules by precluding the stable ex- 
istence of the majority of states in the Hilbert space of the 
system. The physics of decoherence has been thoroughly 
studied during the last few years both from the theoret- 
ical Jl| and also from the experimental point of view [^) . 
As part of these studies it has been recognized that the 
decoherence process has rather unique properties for sys- 
tems whose classical analogues are chaotic ||. In fact, for 
classically chaotic systems decoherence has two very im- 
portant features: First, decoherence restores the validity 
of the correspondence principle that for isolated chaotic 
systems could be violated in an absurdly short time-scale 
(see the Hyperion paradox in Q|). Thus, environment in- 
duced decoherence ensures that classical and quantum 
expectation values agree for a long range of times push- 
ing towards infinity the "breakdown' time which would 
otherwise be present and would depend logarithmically 
on h. This important feature was first conjectured in ||] 
and later studied in more detail in @,||J|. 

The second salient feature of decoherence for classically 
chaotic systems concerns the entropy production rate. In 
fact, in j| it was also conjectured that, when decoher- 
ence is effective, there is a robust range of parameters 
for which the entropy production rate becomes indepen- 
dent of the strength of the coupling to the environment. 
In such regime, the rate is equal to the sum of positive 
Lyapunov exponents (the conjecture was modified in 
where it was noticed that, in general, one should use ap- 
propriately averaged Lyapunov exponents). 

In this letter we present the first conclusive numeri- 
cal evidence supporting the above conjecture for open 
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quantum systems evolving continuously in time (in the 
last few years the status of the conjecture was analyzed 
in part for some open quantum maps |§,§|0). Our re- 
sults show that Lyapunov exponent determines the value 
of the entropy production rate after an initial transient. 
We also clarify other important issues: First we show 
that during the initial transient the entropy production 
rate is proportional to the system-environment coupling 
strength. Second we analyze the transition between the 
two regimes. Our numerical results are consistent with a 
linear dependence of the transition time t c on the entropy 
of the initial state and with a logarithmic dependence on 
the coupling strength. Again, remarkably enough, our 
results validate once more the grossly oversimplified but 
very intuitive picture that was presented in H . 

The system we analyze is a quantum particle mov- 
ing in a quartic double well potential under the ac- 
tion of an harmonic driving force. The Hamiltonian 
is H = p 2 /2 + V(x,t) with the potential V(x,t) = 
—Bx 2 + Cx 4 /2 + Ex cos(wi). This driven system has 
been extensively studied |l0[ and it is well known that 
for a wide range of parameters it has chaotic (mixed) be- 
havior. In our studies we chose several sets of parameters 
in such a way that the stroboscopic phase space portrait 
is like the one shown in Figure 1 where one clearly sees 
the coexistence of islands of stability and a chaotic sea. 
The coupling of this system to an environment is mod- 
eled via the simplest (ohmic, high temperature) quantum 
Brownian motion model. In this case, the evolution of 
the reduced density matrix of the system follows a local 
master equation [fblf , which gives rise to the following 
equation for the Wigner function of the system: 

(_-\ \n% 2n 

W = \H< W) DD 1 ' d {2n+l) Vd {2n+1) W 

n<l V >' 

+ 2 7 9 p ( P W) + Dd 2 pp W (1) 

The physical effects included in this equation are well 
understood. The first term in the right hand side is the 
Poisson bracket that induce classical evolution (Liouville 
flow). The second term carries the quantum corrections. 
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The last two terms bring the influence of the environ- 
ment generating, respectively, friction and diffusion. The 
regime of interest for us is the one where the friction term 
can be neglected (7 — > 0, and energy is approximately 
conserved) while the diffusive term, that produces deco- 
herence is non-negligible (D = 2m7fcsT = constant). 
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FIG. 1. Stroboscopic phase space portrait of the driven 
double well with parameters B = 10, C = 0.5, 
E — l,u> = 5.35. The black ellipses (one in a regular is- 
land and the other in the chaotic sea) show two of the initial 
states we considered (these are pure states and the boundary 
of the ellipses, whose area is 7rln(20)ft ~ 9Ah, corresponds 
to points where the Wigner function decreases to 1/20 of its 
peak value). Throughout our studies we used Ti = 0.1. 

Let us first qualitatively describe the most important 
effects we expect to test in our numerical study. We 
consider initial states with a smooth Gaussian Wigner 
function centered either in the chaotic sea or in a regu- 
lar island as shown in Figure 1 (these two cases lead to 
drastically different results). Let us suppose first that 
the diffusion term in (|I|) is absent (i.e., the state evolves 
according to Schrodinger equation). For a smooth initial 
state the dominant term in (|I|) is the Poisson bracket. 
As Wigner function initially evolves following nonlinear 
classical trajectories it looses its Gaussian shape and de- 
velops tendrils while folding (this happens exponentially 
fast if the initial state is in the chaotic sea). As a con- 
sequence, gradients of W increase and quantum correc- 
tions in (|l|) become important, inducing oscillations in 
the Wigner function. The decoherence producing term 
in (|l|) can be qualitatively understood as being respon- 
sible of two interrelated effects: On the one hand, the 
diffusion term tends to smooth out the Wigner function 
naturally reducing gradients. This washes out the oscilla- 
tions in the Wigner function suppressing quantum inter- 
ference. The time-scale characterizing the dissapearence 
of the fringes can be estimated using previous results [jl3| : 
Fringes with a characteristic wave-vector (along the p— 
axis of phase space) k p decay exponentially with a rate 



given by To = Dk 2 . Noting that a wave-packet spread 
over a distance A x with two coherently interfering pieces 
generate fringes with k p = A x /h one concludes that the 
decoherence rate is = DA 2 /h 2 . This rate, depends 
linearly on the diffusion constant, and is one of the fac- 
tors that contribute to the rate of entropy increase during 
the first phase of the evolution. 

There is a second related aspect of decoherence which 
is drastically different for regular and chaotic systems. 
Thus, appart from suppressing the fringes, the diffusion 
term also tends to spread the regions where the Wigner 
function is possitive, contributing in this way to the en- 
tropy growth. But, as discussed in HJt]], the rate of en- 
tropy production distinguishes regular and chaotic cases. 
For regular states, decoherence should produce entropy 
at a rate which depends on the diffusion constant D. 
However, for chaotic states the rate should become inde- 
pendent of D and should be fixed by the Lyapunov ex- 
ponent. The origin of this /^-independent phase can be 
understood using a simple minded argument (presented 
first in [|| and later discussed in a more elaborate way 
in 0): Chaotic dynamics tends to contract the Wigner 
function along some directions in phase space competing 
against diffusion. These two effects may balance each 
other giving rise to a critical width below which Wigner 
function cannot contract. This local width should be 
approximately a 2 c = 2D/X (being A the local Lyapunov 
exponent). Once this critical size has been reached, the 
contraction stops along the stable direction while the ex- 
pansion continues along the unstable one. Therefore, the 
area covered by the Wigner function grows exponentially 
in time and, as a consequence, entropy grows linearly 
with a rate fixed by the Lyapunov exponents. 

Here, we will present solid numerical evidence support- 
ing the existence of this D-independent phase. For sim- 
plicity, instead of looking at the von Neuman entropy 
T~Cy n = —Tr(p r logp r ) we examine the linear entropy, 
defined as H = — \og(Tr (p 2 )), which is a good measure 
of the degree of mixing of the system and sets a lower 
bound on TLvn- The above argument concerning the 
role of the critical width a c may appear as too simple 
but captures the essential aspects of the dynamical pro- 
cess. Indeed, the master equation can be used to show 
that 

H = 2D((d p W) 2 )/(W 2 ) (2) 

where the bracket denotes an integral over phase space. 
The right hand side of this equation is proportional to the 
mean square wave-number computed with the Fourier 
transform of the Wigner function. This implies that 
the entropy production rate is closely related to the 
phase space structure present in the Wigner distribution. 
Thus, the D-independent phase begins at the time when 
the mean square wave-length along the momentum axis 
scales with diffusion as \[D (as a c does). This behav- 
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ior cancels the diffusion dependence of H which becomes 
entirely determined by the dynamics. 

Appart from analyzing the D-independent phase of 
entropy production we will analyze the nature of the the 
transition time between the diffusion dominated and the 
chaotic regime. This time t c can also be naively esti- 
mated along the lines of the previous argument: The 
time for which the spread of the Wigner function ap- 
proaches the critical one is t c sa A -1 log(a p (0)/a c ). Thus, 
one expects t c to depend logarithmically on the diffusion 
constant and on the initial spread of the Wigner function 
(for Gaussian initial states the spread depends exponen- 
tially on the initial entropy, therefore t c should vary as 
a linear function of the initial entropy) . Our numerical 
work is devoted at testing these intuitive ideas. 

We solved numerically the master equation using two 
different methods (that give similar results). First, we 
computed the Wigner function numerically solving equa- 
tion ([!]) using a third order pseudo-spectral method [ 12) . 
For this, we used a phase space grid of variable size and 
adaptive time-steps (a similar method was used in ||). 
Second we obtained approximate solutions to (|l|) using 
a different strategy: We computed the Floquet states of 
the original (isolated) driven system and obtained the 
master equation in the Floquet basis p]]. This equation 
has a rather simple form when averaged over one driv- 
ing period. In such case, it can be numerically solved by 
usual means (the limiting factor is the number of Floquet 
states that are required). In what follows we present our 
results for the entropy behavior (plots are obtained using 
the first of the two numerical methods, more details and 
other results will be presented elsewhere p"5|). 

The drastic difference between the behavior of the en- 
tropy production rate for regular and chaotic initial con- 
ditions is clearly displayed in Figure 2. Regular initial 
conditions produce entropy at a rate linearly dependent 
on the diffusion coefficient. This is precisely what we 
expect when the entropy production is due to: (i) the 
destruction of interference fringes, (ii) the slow increase 
in the area covered by the possitive Wigner function. The 
oscillations evident in Figure 2 have some distinguishing 
features: they have the frequency of the driving force and 
both their amplitude and their phase are D-independent 
(this is true for the regular initial state and for the initial 
transient of the chaotic state). These oscillations can be 
shown to be related both to the change in orientation of 
the fringes (decoherence is more effective if fringes are 
aligned along the p axis) and to the change in spread 
of the Wigner function in the momentum direction in- 
duced by the dynamics. When the initial state is cen- 
tered in the chaotic sea, the initial transient is followed 
by a regime where the rate 7i becomes independent of 
the value of the diffusion constant (if D is not too small, 
see below). Moreover, the numerical value of the rate 
oscillates around the average local Lyapunov exponent 
computed by averaging over an ensemble of classical tra- 



jectories weighted by the initial Wigner function (In Q| 
the use of a more elaborate averaging scheme was sug- 
gested, we will compare these approaches elsewhere Eqj). 
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FIG. 2. Entropy production rate vs time (in units of the 
driving period). The bold curve is the (time dependent) Lya- 
punov exponent. The linear dependence of the rate on D 
appears in the graph at the top (regular initial state) and 
during the initial transient in the lower plot. In this case (ini- 
tial state in the chaotic sea) the rate becomes independent on 
diffusion and is equal to the Lyapunov exponent (if D the is 
not too small, see text). 

It is remarkable that for long times the entropy pro- 
duction rate is indeed fixed just by the dynamics becom- 
ing independent of D (after all, the entropy production 
is itself a consequence of the coupling to the environ- 
ment but the value of the rate becomes independent of 
it!). The result presented in Figure 2 was shown to be 
robust under changes of initial conditions and other pa- 
rameters characterizing the classical dynamics. There 
are two limitations for the above results to be obtained. 
On the one hand the diffusion constant cannot be too 
strong: In that case the system heats up too fast and 
entropy saturates, making the numerical simulation un- 
reliable. On the other hand, diffusion cannot be too 
small cither: If that is the case decoherence could be- 
come too weak and the interference fringes could persist 
over many oscillations (the minimal value of D required 
for efficient decoherence could be estimated as follows: if 
the Wigner function is coherently spread over a region of 
size A x w 10, we would need a diffusion constant larger 
than D m i n ~ fi 2 1 'A^ ~ 10 -4 for the environment to be 
able to wash out the smallest fringes in one driving pe- 
riod. In fact, Figure 2 shows that when D = 10~ 5 the 
entropy production rate is one order of magnitude smaller 
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than the one corresponding to D = 10 
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FIG. 3. The transition time between the diffusion domi- 
nated regime and the one where the entropy production rate 
is set by the Lyapunov exponent is shown to depend linearly 
on entropy (top) and logarithmically on the diffusion constant 
(bottom). Numerical results were obtained using the param- 
eters: B = 10, C = 0.5, E = 10, w = 6.16, D = 10~ s (top), 
H(0) = (bottom) 

It is interesting to remark that while the simple picture 
presented in ||] is in good qualitative agreement with our 
results, the arguments presented in that paper are too 
simple to include some important effects we found. In 
particular, the oscillatory nature of the rate was com- 
pletely overlooked in However, having said this, it 
is still possible to test some simple results obtained in 
for the transition time between the diffusion domi- 
nated regime and the one where the rate is fixed by the 
Lyapunov exponent. For this purpose we performed two 
related studies. First we analyzed the dependence of the 
transition time t c on the diffusion constant. Of course, 
there is some ambiguity in the definition of the transition 
time because of the oscillatory nature of the rate. Here, 
we defined t c as the time for which the rate reaches some 
value after the initial transient (fortunately the value of 
the rate before and after the transient differ by two or- 
ders of magnitude making this a reasonable definition). 
Second, we analyzed the behavior of the rate as a func- 
tion of the entropy of the initial state (i.e., on the spread 
of the initial Wigner function). Again, we used the same 
definition of t c (for this study we considered a different 
set of parameters than the one used to construct Figure 
1 to assure that the chaotic sea can accommodate ini- 
tial states with linear entropy up to Ti — 4). With this 
results we obtained the plots in Figure 3 where the lin- 



ear dependence of t c on the initial entropy (mixed initial 
states produce entropy at a smaller rate) and its loga- 
rithmic dependence on the diffusion constant are shown. 

Entropy production, a consequence of the entangling 
interactions between the system and its environment, 
plays a central role in studies of decoherence (for exam- 
ple, it is esential to find the "pointer states" of the system 
using the "predictability sieve" ]l6|]). Here, we presented 
results supporting the point of view stating Q that en- 
tropy production rate during decoherence could also be 
used as a diagnostic for quantum chaos. On the other 
hand these results also make evident the fact that, for 
quantum systems that are classically chaotic, the nature 
of the classical limit induced by decoherence is rather pe- 
culiar (and quite different from the one corresponding to 
regular systems). Indeed, this limit exhibits an unavoid- 
able source of unpredictability, being the rate at which 
information is lost into the environment entirely fixed by 
the chaotic nature of the Hamiltonian of the system. 

This work was supported by grants from Anpcyt 
(PICT 01014), Ubacyt and Conicet. 



[1] see for example W. H. Zurek, Physics Today 44, (1991) N 
10, 36; Giulini et al Decoherence and the appearance of the 
classical world in quantum theory (Springer Verlag, Berlin, 
1996). 

[2] J. M. Brune et al, Phys. Rev. Lett. 77, 4887 (1996), see 

also Phys Today51, 7, 36 (1998). 
[3] W. H. Zurek and J. P. Paz, Phys. Rev. Lett. 72 (1994) 

2508. 

[4] W. H. Zurek and J. P. Paz, Physica 83D (1995) 300. 

[5] S. Habib, K. Shizume and W. H. Zurek, Phys. Rev. Lett. 

80 (1998) 4361. 
[6] J. Gong and P. Brumer, Phys. Rev. E60 (1999) 1643. 
[7] A. K. Pattanayak, Phys. Rev. Lett. 83 (1999) 4526. 
[8] K. Shiokawa and B. L. Hu, Phys. Rev E52 (1995) 2497. 
[9] P. A. Miller and S. Sarkar, Phys. Rev. E58 (1998) 4217; 

60 (1999) 1542. 

[10] W. A. Lin and L. E. Ballentine, Phys. Rev. Lett. 65 (1990) 

2927; Phys. Rev.A45 (1992) 3637. 
[11] T. Dittrich, B. Oelschlagel and P. Hanggi, Europh. Lett. 

22 (1993) 5. 

[12] M. Feit, J. Fleck Jr. and A. Steiger, J. Comp. Phys. 47 
(1992) 412. 

[13] J. P. Paz, S. Habib and W. H. Zurek, Phys. Rev. D47 
(1992) 488. 

[14] A. O. Caldeira and A. J. Leggett, Physica 121A (1983) 
587; B. L. Hu, J. P. Paz and Y. Zhang Phys. Rev. D45 
(1992) 2843. 

[15] D. Monteoliva and J. P. Paz, in preparation 

[16] W. H. Zurek, S. Habib and J. P. Paz, Phys. Rev. Lett. 70 

(1993 1187; J. P. Paz and W. H. Zurek, Phys. Rev. Lett. 

82 (1999) 5181. 



4 



